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An analytical approximation of the luminosity distance 
in flat cosmologies with a cosmological constant 

Masaru Adachi and Masumi Kasai*) 

Graduate School of Science and Technology, Hirosaki University, 
Hirosaki 036-8561, Japan 

We present an analytical approximation formula for the luminosity distance in spatially 
flat cosmologies with dust and a cosmological constant. Apart from the overall factor, 
the effect of non-zero cosmological constant in our formula is written simply in terms of 
a rational function. We also show the approximate formulae for the Dyer-Roeder distance 
(empty beam case) and the generalized angular diameter distance from redshift z\ to Z2, 
which are particularly useful in analyzing the gravitational lens effects. Our formulae are 
widely applicable over the range of the density parameter and the redshift with sufficiently 
small uncertainties. In particular, in the range of density parameter 0.3 < J? m < 1 and 
redshift 0.03 < z < 1000, the relative error for the luminosity distance by our formula is 
always smaller than that of the recent work by Wickramasinghe and Ukwatta (2010). Hence, 
we hope that our formulae will be an efficient and useful tool for exploring various problems 
in observational cosmology. 

§1. Introduction 

Current cosmological observations indicate that the universe is spatially flat 
with a cosmological constant. In such a universe, calculations of the luminosity 
distance and the angular diameter distance require numerical integrations and elliptic 
functions^ In order to simplify the repeated numerical integrations, PenP has 
developed an efficient fitting formula. Recently, Wickramasinghe and UkwattsP 
have shown another analytical method, which runs faster than that of Pen® and has 
smaller error variations with respect to redshift z. (See also Ref. H]).) 

In this paper, we present yet another analytical approximation to calculate the 
luminosity distance as follows: 

d L (z, n m ) = U{x(0, fl m )) - -= =<P(x{z, Q m ))\ , (1-1) 

_ 1 + 1.320a + 0.4415x 2 + 0.02656a; 3 
^ ~~ 1 + 1.392x + 0.5121a; 2 + 0.03944a; 3 ' ^ ' ' 

where c is the speed of light, Hq is the Hubble constant, Q m is the density parameter 
of dust matter, related to the density parameter of vacuum energy Qa by i? m + l?yi = 
1, and 

x(z a) = 1 a +v (i ' 3) 

Apart from the overall factor l/\/Q m , the effect of non-zero cosmological constant 
in our distance formula is written simply in terms of a rational function <&(x). 
The function <P{x) has the following properties: 
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1. # = 1 for x = 0. 

2. ^ is a monotonically decreasing function of x, d<P/dx < 0, and <P — >■ 0.6735 for 

x — >■ oo. 

3. (P is a monotonically increasing function of i? m , d<P/df2 m > 0, and ^ — >■ 1 for 
J? m -»• 1. 

4. is a monotonically increasing function of z, d<fr/dz > 0, and <P — > 1 for z — > oo. 

5. 0.6735 < 0(x(O, n m )) < 0(x(z, tt m )) < 1 for < z, < f2 m < 1. 

Note that our approximate formula is explicitly shown to be exact when Q m = 1: 

rfL(M) = ji (1+2) { 1 __y. (1 . 4) 

§2. Approximation 



The luminosity distance in flat cosmologies with a cosmological constant is given 

by 

dL{z,f} m ) = -jjr(l + z) I J"" =. (2-1) 



— z y/f2 m a + (1 - l? m )a 4 
We define ^ 



/o V^a' + (1 - tt m )a' 4 
The power series expansion of F with respect to a around a = yields 



(2-2) 



f =^H*+^ 2 -iI* 3+ ''')' (2 - 3) 

where 

x = l —^a\ (24) 

After expanding F up to 0(x 6 ), we can obtain the Pade approximant to the following 
order: 

r 2 + bix + 6 2 ^ 2 + b 3 x 3 
F = Va— 3— 3, (2-5 

1 + C\X + C2X Z + C^X* 

where the numerical constants are determined as follows: 
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b 2 = 
b 3 = 



C2 



4222975319 


(2-6) 




1599088274 ' 


1138125153117 


(2-7) 


1288865148844' 


7433983569773 


(2-8) 


139933930445920 ' 


635916643 


(2-9) 




L ~ 456882364 ' 


14505955555 


(2-10) 


28326706568 ' 
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44686179629 
1133068262720' 



(2-11) 



Setting a = 1/(1 + z), we finally obtain Eqs. (fl^ - ffbB"]) . 

In order to compare our method with that of Ref . ED , we calculate the following 
relative error 



AE 



I jappr 



d num| 



d num 



100 (per cent), 



(2-12) 



Wickramasinghe & Ukwatta 2010 
This paper 



where d^ pr and df^ m represent the values of luminosity distances calculated by using 
approximate formula and numerical method, respectively. 

Fig. [1] shows a comparison of AE 
for both analytical methods for f2 m = 
0.3. It is shown that our method has a 
smaller relative error for redshift range 
0.03 < z < 100. Although the relative 
error for our method is slightly worse in 
the range 0.01 < z < 0.03, it is still less 
than 0.3 per cent. 

Since our method is based on the 
Taylor expansion and the Pade approx- 
redshift imant with respect to x = a 3 (l — 

Q m )/Q m , it is evident that the error 
in our method decreases monotonically 
with increasing redshift z (or increasing 
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Fig. 



1. A comparison of the percentage rela- 
tive error AE for two analytical methods 
as a function of the redshift for Q m = 0.3. 



f2 m ), i.e., with decreasing x. 



Table I. The percentage relative error AE(%) for the luminosity distance by our formula Eq. 1 ■ 1 p 





z = 0.03 


z = 0.1 


z = l 


2 = 10 


z = 1000 


0.2 


1.87% 


1.38% 


0.25% 


0.08% 


0.06% 


0.3 


0.26% 


0.18% 


0.02% 


0.01% 


< 0.01% 


0.4 


0.03% 


0.01% 


< 0.01% 


<0.01% 


< 0.01% 


1 


















Table [I] shows the relative percentage error AE in our method. It is apparent 
that our approximate formula has sufficiently small uncertainties in the wide range 
of parameters. Only one exception is the nearby region (say, z < 0.1) in the low 
density (say, J? m < 0.2) universe, where the relative error AE exceeds 1 per cent. 
In such a nearby region {z <C 1), however, we may alternatively use the power series 
expansion around z = with sufficient accuracy. 

In particular, in the range of density parameter 0.3 < J? m < 1 and redshift 
0.03 < z < 1000, the relative error for the luminosity distance in our formula is 
always smaller than that of the recent work by Wickramasinghe and UkwattaP* 



4 



M. Adachi and M. Kasai 



§3. The empty beam case 

The distance formula which takes the effect of clumpy distribution of matter 
into account has proposed in Refs. [5]) ,[6j) , and later by Dyer and Roeder^ in more 
general form, which is now known as the Dyer-Roeder distance. The Dyer-Roeder 
luminosity distance for empty beam case in flat cosmologies is 

D L (z,^) = ^-(l + zf f a ] da (3-1) 

Our analytic formula for the empty beam is 



D L l : . <> m ! - l4r ^ + rS I *(*(0. Am)) 1 —9{x[z, Q m )) \ , (3-2) 

(1+Z)2 J 



5 Hn v fi„ 



_ 1 + 1.256a + 0.3804s 2 + 0.0164a 3 
^ ~ 1 + 1.483s + 0.6072s 2 + 0.0587x 3 ' ^ ' ' 



where s is defined in Eq. (|1 -3)) . The formula can be obtained in the same way as that 
in the previous section. The rational function \P(x) has the following properties: 

1. \P = 1 for x = 0. 

2. & is a monotonically decreasing function of x, d\P/dx < 0, and & — > 0.2792 for 

x — > CO. 

3. ^ is a monotonically increasing function of i? m , d\P/df2 m > 0, and IF — >• 1 for 
^ m -> 1. 

4. •f' is a monotonically increasing function of z, d\P/dz > 0, and 1 for 

Z — >• CO. 

5. 0.2792 < «f(s(0, ft m )) < !P(a?(*, ft m )) < 1 for < z, < l? m < 1. 
Again, our approximate formula Eq. (|3-2|) is exact when Q m = 1: 

D h (z, 1) = H ° (i + z) 2 (l - -^-t) • (3-4) 
o^o [ (l + z)2 J 

Note that no consideration on the Dyer-Roeder distance has been taken in the pre- 
vious papers, e.g., Refs. [2]) ,[3]). 

Table [II] shows the relative error of 

T U1 TT rrii , i Anfws the our formula Eq. (13-20 . The relative 
Table II. The percentage relative error AE(%) M 11 ' 

for the Dyer-Roeder distance (empty beam error for 0.3 < J? m < 1 is always less 

than 0.15 per cent in the range 0.03 < 
z < 10. For i? m = 0.2, the accuracy 
gets slightly worse, but the error is still 
less than 1 per cent in the same redshift 
range. 

We omitted the error calculations 
for z > 10 for the following reasons. 
First, the errors are sufficiently small in 
those regions, and second, the Dyer-Roeder description is relevant only in the regions 
where the clumpy distribution of matter becomes important. 



case) 


by our 


formula Eq. 


(EH) 




iT2 m z 


= 0.03 


z = 0.1 


z = 1 


z = 10 


0.2 


0.90% 


0.68% 


0.21% 


0.15% 


0.3 


0.14% 


0.10% 


0.03% 


0.02% 


0.4 


0.02% 


0.01% 


< 0.01% 


< 0.01% 


1 
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§4. The generalized angular diameter distance 



Here we present the analytic formulae for the generalized angular diameter dis- 
tance from redshift z = Z\ to z = Z2, which is frequently required in analyzing the 
gravitational lens effects. Here we only consider the case of flat cosmologies. For the 
distance formulae in non-flat cases, see, e.g., Ref. [8]). For the standard (filled beam) 
case, the angular diameter distance from redshift z = Z\ to z = Z2 is 

Our analytical approximation is simply 
d A (zi, z 2 ) = 7==^ { J— ${x{z u Q m )) - J— $(x(z 2 , f2 m )) 

-^0 V^m(l + Z2) y VI + z l V 1 + z 2 

(4-2) 

where <&{x) is defined in Eq. (|1-2|) . The relative error of our formula Eq. (|4-2p for 
f2 m = 0.3 is less than 0.02 per cent in the range < Z\ < 1 for z 2 = 1, and less than 
0.01 per cent in the range < z\ < 3 for z 2 = 3. 
For the empty beam case, 

j — / . c . f 1+21 a 2 da .. „. 

Da (zuZ2 ) = jr I + / ==, 4-3 

and our approximate formula is 



5#0 V^m [(1 + ^)2 (1 + ^)5 J 

(4-4) 

where ^(x) is already defined in Eq. (|3-3|) . 

The reciprocity theorem holds between the angular diameter and luminosity 
distances as follows: di,(z) = (1 + z) 2 d A {Q, z), and D^(z) = (1 + z) 2 Da(0, z). 



§5. Summary 



We have presented a simple analytical approximation formula for the luminosity 
distance in flat cosmologies with a cosmological constant. We have also shown the 
approximate formulae for the Dyer-Roeder distance and the generalized angular 
diameter distance from redshift z = z\ to z = z 2 , which are particularly useful 
in analyzing the gravitational lens effects. Apart from the overall factor l/\/f2 m , 
the effects of non-zero cosmological constant in our distance formulae are written 
simply in terms of the rational functions <£(x) for "filled beam case" and H/(x) for 
"empty beam case". Both are monotonically decreasing functions with respect to x, 
and increasing ones with respect to redshift z and the density parameter fi m . 

Our formulae are widely applicable over the range of the density parameter 
and the redshift with sufficiently small uncertainties. In particular, in the range 
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0.3 < i? m < 1 and 0.03 < z < 1000, the relative error for the luminosity distance by 
our formula is always smaller than that of the recent work by Wickramasinghe and 
Ukwatta.® Hence, we hope that it will be an efficient and useful tool for exploring 
various problems in observational cosmology, such as the statistical grativational 
lensing, the cosmological parameter fitting in the magnitude-redshift relation of the 
super novae, and so on. 

Appendix A 

Approximations in a small redshift region 



In the region of small redshift {z <C 1), the standard power series expansions 
can safely be used. The power series expansions of d^(z, i? m ) (filled beam case) and 
Dl(z, i? m ) (empty beam case) around z = are 

d L (z, f2 m ) = ^ L + (l - z* + 9 ^~ W Q m z 3 + •••}, (Al) 

D L {z, D m ) = 4-\ z +( 1 - 7^m J z 2 + r~Q m - 1 ) n m z 3 + •••}. (A-2) 



Ho I V 4 
Just for reference, their Pade approximants are 

c (12f? m -16)z + (9^ + 4^ m -16)z 2 
U ' m) ~ H (12fl m - 16) + (18/4 " 20^ m )z ' 1 j 

c (12^ m -16)z + (9^ + 8^ m -16)z 2 
Ll ' mj " i?o (12fi m - 16) + (18^ - I6fl m )z ' 1 j 



Table III. The percentage relative error AE(%) of the power series expansion Eq. (JA-1 



iT2 m 


z = 0.01 


2 = 0.1 


z = 0.2 


z = 0.5 


2 = 1 


0.2 


< 0.01% 


< 0.01% 


0.02% 


0.10% 


0.12% 


0.5 


< 0.01% 


0.01% 


0.09% 


1.25% 


8.04% 


0.7 


< 0.01% 


0.02% 


0.15% 


1.80% 


10.56% 


1.0 


< 0.01% 


< 0.01% 


0.05% 


0.66% 


3.98% 



Table IV. The percentage relative error AE(%) of the Pade approximant Eq. (|A-3fl 



£2m 


z = 0.01 


z = 0.1 


2 = 0.2 


2 = 0.5 


2=1 


0.2 


< 0.01% 


< 0.01% 


0.05% 


0.51% 


2.30% 


0.5 


< 0.01% 


< 0.01% 


0.03% 


0.24% 


0.72% 


0.7 


< 0.01% 


< 0.01% 


< 0.01% 


0.05% 


0.36% 


1.0 


< 0.01% 


< 0.01% 


< 0.01% 


0.09% 


0.42% 



The relative error of Eqs. (|ATj) and (|A-3p are listed in Tables ITO1 and HVl The 
maximal relative error in the power series expansion AE is 10.6 per cent at z = 1 in 
the interval < z < 1 and 0.2 < f2 m < 1.0, not 37 per cent which was claimed by 
PenP The Pade approximant Eq. (|A-3j) shows in many cases better accuracy than 
the power series expansion Eq. (|A-ip , The relative error of the Pade approximant 
does not exceed 3 per cent even at redshift z = 1 in the rage 0.2 < i? m < 1.0. 
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Appendix B 

An analytical approximation of the growth function 

There is one more thing. Here we present another approximate formula which 
we hope to be an efficient and useful tool in observational cosmology. Heath^ has 
shown that the growth function in a dust cosmology can be written as 



Di(a) oc H(a) / / — (B-l) 



da' 

lo (a'H(a')Y 
D 2 <xH(a), (B-2) 

H{a) = v^ma -3 + (1 — J7 m — Q A )a~ 2 + Q A . (B-3) 

Although a compact expression using the incomplete beta function has been 
shown in Ref. I1U|) . so far, no analytic solution of D\{a) has been presented for 
Qa 7^ 0. Here we restrict ourselves to the case Q m + Qa = 1, and present an 
approximate formula in a simple algebraic form. 

We adopt a normalization for D\{a) as 

Dl (a) = ^H(a) [ a -^— 3 . (B-4) 
z Jo {a'H(a')) 

Then, in a flat cosmology Q m + Q A = 1, our formula is 

, , 1 + 1.175s + 0.3064s 2 + 0.005355s 3 ,„ . 

D\(a) = ay 1 + x k 3 — , B-5) 

u 1 1 + 1.857s + 1.021x 2 + 0.1530s 3 ' v ; 



where, 



^ ^ m 3 fr> a\ 

x = — — a . (B"6) 



Note that our approximate formula is exact when J7 m = 1. 

The well-known approximation formula for the growth function in Ref. [TT|) . 
which was adopted from Ref. I12|) , is 

D? = (B-7) 

for f2 m + Qa = 1. A comparison of the relative error AE at a = 1 is shown in 
Table IVl for Q m + Qa = 1- Our formula has generally smaller relative error in the 
range 0.2 < Q m < 1. 



Table V. A comparison of the relative error AE(%) of Eq. (|B-5|) and Eq. (|B-7|) at a = l(z = 0). 



iT2 m 




0.2 


0.3 


0.5 


0.9 


AE of Eq. 
AE of Eq. 


CEP 

(lB~7ll 


0.19% 
0.54% 


< 0.01% 
0.134% 


< 0.01% 
0.057% 


< 0.01% 
0.019% 
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